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Abstract 

Background: Color and pH of meat measured 24 h post mortem are common selection objectives in pig breeding 
programs. Several amino acid substitutions in PRKAG3 have been associated with various meat quality traits. In our 
previous study ASGA0070625, a SNP next to PRKAG3, had the most significant association with meat quality traits in 
the Finnish Yorkshire. However, the known amino acid substitutions, including 1199V, did not show any association. 
The aims of this study were to characterize further variation in PRKAG3 and its promoter region, and to test the 
association between these variants and the pH and color of pork meat. 

Results: The data comprised of 220 Finnish Landrace and 230 Finnish Yorkshire artificial insemination boars with 
progeny information. We sequenced the coding and promoter region of PRKAG3 in these and in three additional 
wild boars. Genotypes from our previous genome-wide scans were also included in the data. Association between 
SNPs or haplotypes and meat quality traits (deregressed estimates of breeding values from Finnish national 
breeding value estimation for pH, color lightness and redness measured from loin or ham) was tested using a linear 
regression model. Sequencing revealed several novel amino acid substitutions in PRKAG3, including K24E, 141V, 
K131R, and P134L. Linkage disequilibrium was strong among the novel variants, SNPs in the promoter region and 
ASGA0070625, especially for the Yorkshire. The strongest associations were observed between ASGA0070625 and 
the SNPs in the promoter region and pH measured from loin in the Yorkshire and between 1199V and pH measured 
from ham in the Landrace. In contrast, ASGA0070625 was not significantly associated with meat quality traits in the 
Landrace and 1199V not in the Yorkshire. Haplotype analysis showed a significant association between a haplotype 
consisting of 799/ and 24E alleles (or g-157C or g-58A alleles in the promoter region) and pH measured from loin 
and ham in both breeds (P-values varied from 1.72 x 10~ 4 to 1.80 x 10~ 8 ). 

Conclusions: We conclude that haplotype g-157C - g.-58A - 24E - 1991 in PRKAG3 has a positive effect on meat 
quality in pigs. Our results are readily applicable for marker-assisted selection in pigs. 
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Background 

Meat quality characteristics such as water-holding cap- 
acity, tenderness, intramuscular fat, and taste are im- 
portant for the meat industry as well as for consumers 
[1-3]. Since direct measurement of some of these traits 
is considered laborious in practice, many breeding pro- 
grams only use correlated traits like pH and color of 
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meat for selection. In Finland, pH and color of loin have 
been among the pig breeding objectives since 1983 and 
meat quality of ham since 2000 [1]. 

Several important genes are known to have a major effect 
on meat quality and carcass composition traits in pigs, 
including RYR1 (ryanodine receptor 1) on chromosome 6 
[4]; PRKAG3 (protein kinase, AMP-activated, gamma 3 
non-catalytic subunit) on chromosome 15 [5,6]; IGF2 (insu- 
lin-like growth factor 2) on chromosome 2 [7,8]; CAST 
(calpastatin) on chromosome 2 [9]; andMC4i? (melanocor- 
tin 4 receptor) on chromosome 1 [10]. A previous genome- 
wide association study in the Finnish Yorkshire population 
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revealed a strong association between pH measured from 
loin and a chromosomal region around PRKAG3 [11]. 
However, none of the reported amino acid substitutions 
T30N, G52S, L53P, I199V, or R200Q were as strongly asso- 
ciated with pH as SNP ASGA0070625 (RS80816788, pos- 
ition 133,677,385 Sus Scrofa build 10.2), which is located 
near the PRKAG3 gene. A further study with more markers 
on the PRKAG3 region was therefore needed. Additionally, 
the Finnish Landrace breed could be used as a validation 
population for the initial findings. 

This article reports novel variations within the 
PRKAG3 gene in Finnish Yorkshire and Landrace pigs, 
shows the linkage disequilibrium (LD) structure around 
PRKAG3, and presents the results from our association 
analysis between SNPs or haplotypes and meat quality 
traits. We show that to achieve a positive effect on meat 
quality traits, particularly in pH measured post mortem 
from loin and ham, the animal has to carry both the 
1991 allele and the combination g.-157C - g.-58A - 24E. 
This haplotype was significant for meat quality in both 
breeds. 

Results 

SNP discovery 

Sequencing of the exons and promoter region of 
PRKAG3 revealed four novel amino acid substitutions: 
K24E, 141 V, K131R, and P134L (Table 1). The four SNPs 
upstream of the PRKAG3 transcription initiation site re- 
ported by Ryan et al. [12] were also detected in the stud- 
ied Finnish Yorkshire and Landrace populations. Two 
additional synonymous substitutions were observed at 
amino acids 193 and 194, but these were very rare in 



both populations (frequency less than 3%) and were thus 
excluded from later analysis. Minor allele frequencies of 
the novel amino acid substitutions varied from 0.17 to 
0.20 and from 0.15 to 0.20 in the Yorkshire and 
Landrace, respectively. Alleles 53P and 200R were fixed 
in both populations. Allele 30T was also very rare in the 
Yorkshire, with a frequency of only 1% compared to 13% 
in the Landrace, whereas alleles S2S and 1991 were more 
common in the Yorkshire than in the Landrace. 

Functional analysis of the amino acid substitutions using 
SIFT [13] (Sorting Intolerant From Tolerant, http://sift.bii. 
a-star.edu.sg/) underlined three substitutions, K24E, I41V 
and L53P, as damaging (Table 2). These effects appeared to 
be transcript specific indicating differences between effects 
of substitutions on protein isoforms. For substitutions I41V 
and L53P the SIFT scores were similar for all protein iso- 
forms, but a clear difference for K24E was identified 
(Table 2). 

Linkage disequilibrium 

Figure 1 shows the linkage disequilibrium (LD) within 
and around PRKAG3. Overall, LD was stronger and ex- 
tended over a longer range in the Finnish Yorkshire than 
Landrace. In the Yorkshire, g.-311A > G, g.-221G > A, g.- 
157C>G, g.-58A>G, K24E, 141 V, K131R, and P134L 
were in complete LD with ASGA0070625, the SNP that 
showed the strongest association with pH measured 
from loin in our previous whole-genome analysis [11]. 
In the Finnish Landrace, ASGA0070625 was in complete 
LD with g.-157C>G, g.-58A>G, and K24E. Interest- 
ingly, the known amino acid substitutions T30N, G53S, 
and I199V were in very weak LD with ASGA0070625 



Table 1 Positions, alleles and minor allele frequencies of the identified SNPs in PRKAG3 



Exon 


Variation 3 


SS number 


RS number 


Position on chromosome 15, bp b 


Alleles" 


Yorkshire 


Landrace 




ASGA0070625 




80816788 


133677385 


A/G 


0.21 


0.28 




g.-31 1 A > G 




196959880 


133800456 


G/A 


0.19 


0.13 




g.-221G > A 




196952335 


133800546 


A/G 


0.19 


0.13 




g.-157C > G 




196956394 


133800610 


G/C 


0.20 


0.26 




g.-58A > G 




196959698 


133800709 


G/A 


0.20 


0.26 


2 


K24E 


947848631 




133801207 


A/G 


0.17 


0.20 


3 


T30N 




328566929 


133802071 


A/C 


0.01 


0.13 


3 


141V 


947848632 




133802103 


A/G 


0.20 


0.16 


3 


G52S 




343733804 


133802136 


A/G 


0.43 


0.12 


4 


L53P 




337462352 


133802499 


T/C 


0 


0 


4 


K131R 


947848633 




133802733 


G/A 


0.19 


0.15 


4 


P134L 


947848634 




133802742 


T/C 


0.19 


0.16 


5 


1199V 






133803828 


A/G 


0.45 


0.30 


5 


R200Q 






133803829 


A/G 


0 


0 



a The novel amino acid substitutions are given in italics with corresponding dbSNP submission numbers (http://www.ncbi.nlm.nih.gov/SNP/). 
b Based on Sus Scrofa build 10.2. 
c The minor allele is given first. 
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Table 2 Effect of the identified SNPs in PRKAG3 on protein sequence 



Reference 


Position on chromosome 15, bp 


Amino acid co-ordinate 


Transcript 


SIFT 3 


K24E 


133801207 


47 


ENSSSCT00000017641 


0.16 


K24E 


133801207 


24 


ENSSSCT00000033825 


0.01 


T30N 


133802071 


80 


ENSSSCT00000017641 


0.44 


T30N 


133802071 


57 


ENSSSCT00000033825 


0.20 


T30N 


133802071 


30 


ENSSSCT00000036402 


0.32 


141V 


133802103 


91 


ENSSSCT00000017641 


0.05 


141V 


133802103 


68 


ENSSSCT00000033825 


0.07 


141V 


133802103 


41 


ENSSSCT00000036402 


0.03 


G52S 


133802136 


102 


ENSSSCT00000017641 


0.11 


G52S 


133802136 


79 


ENSSSCT00000033825 


0.13 


G52S 


133802136 


52 


ENSSSCT00000036402 


0.13 


L53P 


133802499 


103 


ENSSSCT00000017641 


0.04 


L53P 


133802499 


80 


ENSSSCT00000033825 


0.05 


L53P 


133802499 


53 


ENSSSCT00000036402 


0.07 


K131R 


133802733 


181 


ENSSSCT00000017641 


1.00 


K131R 


133802733 


131 


ENSSSCT00000036402 


1.00 


K131R 


133802733 


158 


ENSSSCT00000033825 


1.00 


P134L 


133802742 


184 


ENSSSCT00000017641 


0.58 


P134L 


133802742 


134 


ENSSSCT00000036402 


0.58 


P134L 


133802742 


161 


ENSSSCT00000033825 


0.58 


1199V 


133803828 


249 


ENSSSCT00000017641 


1.00 


1199V 


133803828 


199 


ENSSSCT00000036402 


1.00 


1199V 


133803828 


226 


ENSSSCT00000033825 


1.00 



a SIFT score <0.05 is considered as damaging. 



and, hence, also with the novel amino acid substitutions 
and SNPs in the promoter region of PRKAG3 in both 
breeds. 

Haplotypes 

In total, 12 different haplotypes were observed in ei- 
ther of the breeds (Table 3). One of these, HAP3, is 
actually the same as the wild haplotype identified in 
the sequences from the European wild boars used in 
this study (based on three wild boar samples). HAP3 
was the most common haplotype in the Finnish 
Landrace breed (frequency 0.34), but very rare in the 
Finnish Yorkshire (frequency 0.01). HAP1, a haplotype 
similar to the wild haplotype except for a point muta- 
tion of G to A at G52S, was the most frequent in the 
Yorkshire with a frequency of 0.42. Again there was a 
notable difference in frequencies between the two 
breeds (frequency only 0.12 in the Landrace). HAP2, 
another haplotype similar to the wild haplotype ex- 
cept for a point mutation of G to A at 1 199V, was 
common in both breeds, with frequencies of 0.36 and 
0.25 for the Yorkshire and Landrace, respectively. The 
other haplotypes (HAP4 to HAP 12) jointly accounted 



for 22% to 29% of the haplotypes identified in the 
Yorkshire and Landrace, respectively. Some of the 
rare haplotypes may be spurious due to genotyping or 
haplotyping errors. 

SNP association 

Based on the genotypes from our previous whole- 
genome scans [11,14], the most significant SNPs 
affecting meat quality traits were in the region from 
120 Mb to 140 Mb on chromosome 15 in the 
proximity of PRKAG3 (Figure 2). The smallest 
P-value was observed for association between six 
SNPs (ASGA0070634, ASGA0070625, MARC0083357, 
DBUN0002708, MARC0039273, and DIAS0002965) 
and pH measured from loin in the Finnish Yorkshire. 
These six SNPs were in complete LD in the region 
from 133.64 Mb to 134.01 Mb. Their estimated allele 
substitution effect for the Yorkshire was -0.059 ± 
0.008 (P-value = 7.28 x 10" 13 ), corresponding to 1.3 SD 
of the estimated polygenic effect. Minor alleles A, A, 
C, A, A, and G of ASGA0070634, ASGA0070625, 
MARC0083357, DBUN0002708, MARC0039273, and 
DIAS0002965, respectively, decreased pH measured 
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Figure 1 Linkage disequilibrium (LD) expressed as r 2 values (black color indicated complete LD) within and around PRKAG3 for the 
Finnish Yorkshire (A) and Landrace (B). 



from loin. None of these six SNPs showed statistically 
significant association with any of the measured traits 
in the Finnish Landrace (see Figure 2 for pH 
measured from loin, as an example). Statistical 
significance was claimed if the P-value was below 
2.0 x 1(T 6 . This stringent limit was set because the 
observations used in this study were the same as 
those in our previous genome- wide scans [11,14]. 

Because of the strong LD with ASGA0070625, the 
novel amino acid substitutions identified in this study 
and SNPs in the promoter region of PRKAG3 gave 
very small P-values for pH measured from loin for 



the Finnish Yorkshire (Table 4). The differences in P- 
values between ASGA0070625 and the SNPs which 
were in complete LD with ASGA0070625 were due to 
slight differences in the number of genotypes available 
for each SNP. No significant association was detected 
between I199V and the measured meat quality traits 
for the Yorkshire. For the Landrace, in contrast, the 
strongest association was observed between 1 199V 
and pH measured from ham (P-value = 6.44 x 10~ 7 
with an allele substitution effect of 0.030 ± 0.006). The 
association between I199V and pH measured from 
loin was also strong in the Landrace (Table 4). 
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Table 3 Haplotypes in PRKAG3 and their frequencies 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


1 1 


12 


Wild 


ASGA0070625 


G 


G 


G 


G 


A 


A 


A 


A 


A 


A 


A 


A 




g -31 1 A > G 


A 


A 


A 


A 


A 


A 


A 


G 


G 


G 


G 


G 


A 


g.-221G > A 


G 


G 


G 


G 


G 


G 


G 


A 


A 


A 


A 


A 


G 


g.-157C>G 


C 


C 


C 


C 


G 


G 


G 


G 


G 


G 


G 


G 


C 


g -58A > G 


A 


A 


A 


A 


G 


G 


G 


G 


G 


G 


G 


G 


A 


K24E 


G 


G 


G 


G 


A 


G 


G 


A 


A 


A 


A 


A 


G 


I dUN 


r 




c 


L 


A 

A 


A 

A 


A 

A 


L 


r 


L 


L 


A 

A 


L 


141V 


G 


G 


G 


G 


G 


G 


G 


A 


A 


A 


G 


G 


G 


G52S 


A 


G 


G 


A 


G 


G 


G 


G 


G 


G 


G 


G 


G 


K131R 


A 


A 


A 


A 


A 


A 


A 


G 


G 


G 


A 


A 


A 


P134L 


C 


C 


C 


C 


C 


C 


C 


T 


T 


T 


C 


C 


C 


1199V 


G 


A 


G 


A 


G 


A 


G 


A 


G 


G 


G 


G 


G 


Yorkshire 


0.418 


0.363 


0.007 


0.007 


0.002 


0.002 


0.007 


0.080 


0.087 


0.027 


0 


0 




Landrace 


0.122 


0.249 


0.341 


0 


0.124 


0 


0 


0.041 


0.116 


0 


0.002 


0.004 





boar 



Haplotype association 

Four haplotypes were identified having allele G of 
ASGA0070625 (Table 3). When each of these was tested 
against all other haplotypes in the Yorkshire data, the 
only statistically significant association (P-value = 6.35 x 
10" 8 , Table 5) was observed between HAP2 (with 1991 
or allele A) and pH measured from loin. The substitu- 
tion effect of this haplotype was 0.039 ± 0.007. HAP2 
was also the only haplotype that showed a significant as- 
sociation in the Finnish Landrace with both pH mea- 
surements (Table 5). The haplotype substitution effects 




120.00 122.00 124.00 126.00 128.00 130.00 132.00 134.00 136.00 138.00 140.00 



Figure 2 P-values (-loglO) of the SNPs for pH measured from 

loin. P-values for the Finnish Yorkshire are marked with black 

diamonds and for the Finnish Landrace with green diamonds. 
\ ) 



for pH measured from loin and ham were 0.031 ± 0.007 
and 0.039 ± 0.007, respectively. Haplotypes carrying al- 
lele A of ASGA0070625 (HAP5-HAP12) had no signifi- 
cant association with any of the tested meat quality 
traits. Haplotypes with only one or two observations 
(HAP6, HAP11, and HAP12) were excluded from the 
analysis. 

Discussion 

Meat color and pH are traits which are commonly in- 
cluded in pig breeding programs to improve the techno- 
logical properties of pork and to increase consumer 
gratification. Several studies indicate that PRKAG3 is 
one of the key genes causing variation in pork meat pH, 
L* (lightness of color), and drip loss between animals 
[5,6,15-18]. The protein encoded by PRKAG3 is the skel- 
etal muscle cell- specific regulatory subunit gamma 3 of 
AMP-activated protein kinase (AMPK). AMPK is an en- 
ergy sensor which, when activated in response to cellular 
metabolic stresses, directly phosphorylates and inacti- 
vates the key enzymes involved in regulating de novo 
biosynthesis of fatty acid and cholesterol. The best 
known mutation in PRKAG3 is 200Q, which is found 
only in the Hampshire pig breed. The allele 200Q causes 
a high content of stored glycogen in white skeletal mus- 
cles, leading to low muscle pH 24 h post mortem, poor 
water-holding capacity, and low processing yield [5]. 
Furthermore, 1 199V and T30N are reported to affect pH 
[6,15,17,18], and variations in the promoter region of 
PRKAG3 have been associated with gene expression and 
meat quality [12]. 

In this study, we characterized several novel amino acid 
substitutions within exons 2, 3, and 4. The SNPs character- 
ized in the PRKAG3 promoter region correspond to those 
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Table 4 P-values of the association between SNPs in PRKAG3 and meat quality traits 



Yorkshire Landrace 



Variation 






Loin 












Ham 






Loin 






Ham 




PH 




L* 




a* 




PH 




L* 


a* 


PH 


L* 


a* 


PH 


L* 


a* 


ASGA0070625 


7.27 x 10" 


-13 


5.52x10" 


-6 


4.38x10" 


-5 


3.47 x 10" 


-6 


0.73 


8.85 X10" 3 


0.02 


0.03 


0.03 


0.02 


0.34 


0.79 


g.-311A>G 


6.83x10" 


-10 


4.22 x 10" 


-5 


6.28 X10" 4 


1.44x10" 


-4 


0.74 


0.05 


0.28 


0.36 


0.36 


0.60 


0.93 


0.79 


g -221G > A 


6.83x10" 


-10 


4.22 x 10" 


-5 


6.28x10" 


-4 


1.44x10" 


-4 


0.74 


0.05 


0.33 


0.62 


0.36 


0.64 


0.98 


0.77 


g.-157C>G 


7.94 x 10" 


-1 1 


2.79 x 10" 


-5 


4.29x10" 


-4 


2.68 x 10" 


-5 


0.56 


0.03 


0.02 


0.06 


0.37 


0.05 


0.38 


0.83 


g.-58A > G 


7.94 x 10" 


-1 1 


2.79 x 10" 


-5 


4.29 X10" 4 


2.68x10" 


-5 


0.56 


0.03 


0.03 


0.09 


0.36 


0.06 


0.41 


0.92 


K24E 


5.80 x 10" 


-10 


3.83x10" 


-5 


3.16X 10" 


-3 


7.83 x 10" 


-5 


0.38 


0.03 


0.23 


0.27 


0.16 


0.15 


0.57 


0.54 


T30N 


0.13 




0.39 




0.46 




0.12 




0.61 


0.36 


0.02 


6.51 X10" 3 


0.32 


0.01 


0.16 


0.50 


141V 


1.41 x 10" 


-1 1 


1.43 x 10" 


-5 


9.20x10" 


-5 


1.96x10" 


-5 


0.82 


0.02 


0.45 


0.92 


0.04 


0.68 


0.53 


0.59 


G52S 


0.36 




0.58 




0.21 




0.77 




0.91 


0.81 


0.41 


0.79 


0.84 


0.78 


0.53 


0.64 


K131R 


1.41 x 10" 


-1 1 


1.43x10" 


-5 


9.20 x 10" 


-5 


1.96x10" 


-5 


0.82 


0.02 


0.41 


0.93 


0.08 


0.54 


0.60 


0.62 


P134L 


1.41 x 10" 


-1 1 


1.43x10" 


-5 


9.20 x 10" 


-5 


1.96x10" 


-5 


0.82 


0.02 


0.43 


0.85 


0.04 


0.70 


0.55 


0.67 


1199V 


0.01 




0.19 




0.12 




0.09 




0.93 


0.10 


1.62 X10" 5 


0.09 


0.04 


6.44 x 10" 


" 7 0.17 


0.67 


DBUN0002708 


3.21 x 10" 


-10 


4.15 x 10" 


-5 


3.81 x 10" 4 


3.36x10" 


-4 


0.87 


0.01 


5.37 X10" 3 


0.02 


0.05 


0.07 


0.37 


0.75 



reported by Ryan et al. [12]. Based on the genomic se- 
quence the SNP g.-158C > G reported by Ryan et. al. corre- 
sponds to our g.-157C > G. The genomic location of these 
SNPs is the same. All novel amino acid substitutions 
(K24E, I41V, K131R, and P134L) and characterized SNPs in 
the promoter region were in complete LD in the Finnish 
Yorkshire, based on Haploview analysis. The Finnish Land- 
race breed showed more diversity, with only g.-157C > G, 
g.-58A > G, and K24E in complete LD with each other and 
with ASGA0070625. Strong but not complete LD was also 
reported by Ryan et al. [12] for SNPs in the promoter re- 
gion. No significant LD was found between the promoter 
region SNPs or novel amino acid substitutions and I199V. 
This finding is similar to the observation made by Ryan 
et al. [12] regarding LD between promoter region SNPs 
and I199V. 



A slightly different picture of LD can be drawn from 
the haplotype estimates obtained by the FastPHASE pro- 
gram (Table 3). Haplotypes HAP5, HAP6, and HAP7 
divide the LD pattern into two groups. The first group 
comprises the SNPs in complete LD, namely 
ASGA0070625, g.-157C > G, and g.-58A > G in the 
Finnish Yorkshire, and additionally K24E in Finnish 
Landrace. For Landrace, this is exactly the same result 
as given by Haploview. The second group includes SNPs 
g.-311A>G, g.-221G>A, I41V, K131R, and P134L in 
the Yorkshire, but not in the Landrace. 

There are several explanations for the different LD 
outcomes from Haploview and FastPHASE analyses, 
such as different algorithms and ways of treating missing 
genotypes between the two programs. Genotyping or 
phasing errors may also cause spurious haplotypes. The 



Table 5 P-values of the associations between haplotypes in PRKAG3 and meat quality traits 

Yorkshire Landrace 



Haplotype 3 




Loin 






Ham 






Loin 






Ham 




PH 


L* 


a* 


PH 


L* 


a* 


PH 


L* 


a* 


PH 


L* 


a* 


HAP1 


0.47 


0.90 


0.14 


0.76 


0.69 


0.66 


0.41 


0.89 


0.92 


0.65 


0.61 


0.57 


HAP2 


6.35 X10" 8 


6.00 X10" 4 


0.01 


1 .72 x 1 0" 4 


0.97 


0.01 


8.69 X10" 6 


0.15 


3.99 X10" 3 


1.80x10" 


" 8 0.35 


0.41 


HAP3 


0.15 


0.51 


0.17 


0.40 


0.03 


5.98 X10" 3 


0.04 


0.41 


0.17 


0.02 


0.74 


0.13 


HAP4 


0.42 


0.06 


0.24 


0.88 


0.23 


0.38 














HAP5 














0.02 


0.01 


0.27 


0.02 


0.17 


0.60 


HAP7 


0.73 


0.44 


0.78 


0.16 


0.80 


0.32 














HAP8 


2.25 X10" 4 


3.24 x 10" 3 


0.11 


3.03 X10" 3 


0.82 


0.35 


0.12 


0.05 


0.83 


0.12 


0.08 


0.48 


HAP9 


2.29 X10" 4 


0.05 


4.98 x 10" 4 


0.13 


0.48 


0.08 


0.07 


0.23 


0.12 


0.16 


0.63 


0.94 


HAP10 


0.04 


0.17 


0.50 


0.11 


0.38 


0.42 















a Because of low number of observations, haplotypes HAP6, HAP1 1, and HAP12 were excluded from the analysis. 
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latter view is supported by the fact that haplotypes 
HAP6 and HAP7 are extremely rare in the Yorkshire 
and absent in the Landrace, while haplotypes HAP 11 
and HAP 12 are completely absent in the Yorkshire and 
are carried only by one or two animals in Landrace. 
Additionally, some haplotypes may have been introduced 
into one breed from the other through occasional invol- 
untary crossing of breeds at the farm level. 

Single-SNP analysis yielded controversial results when 
the two breeds were compared. The SNP which was re- 
ported as highly significant for pH measured from loin 
in the Finnish Yorkshire in our previous study [11] was 
significant also in this analysis, given the fact that most 
of the animals (Yorkshire boars) were the same in both 
analyses. However, had the Finnish Landrace been used 
as the validation population for the previous study, the 
significance of ASGA0070625 would not have been re- 
peated and this SNP would have been claimed to be a 
population-specific marker for meat pH. Similarly, 
1 199V was not repeated in the Yorkshire, raising a doubt 
that 1 199V is breed- or population-specific. But when 
haplotypes instead of single SNPs were used in the asso- 
ciation analysis, the results were coherent: the same 
haplotype (HAP2) was significantly associated with pH 
in both breeds. 

The haplotype with both 1991 and 24E alleles (or g-157C 
or g.-58A) was found favorable for pH measured from loin 
and ham in both breeds. This provides strong support for 
the hypothesis that allele 1991 alone does not create a posi- 
tive effect on the pH level in muscle post mortem, but the 
animal has to carry an additional variation either in the pro- 
moter region of PRKAG3 (g.-157C or g.-58A or both) or 
glutamate at amino acid position 24 (or 47 depending on 
the PRKAG3 isoform used for naming). Analysis of the 
SNP effect on the protein function suggested that there 
may be some differences between transcripts. Based on a 
SIFT analysis [13], the K24E mutation showed a significant 
(SIFT score < 0.05) effect on the protein function in the 
ENSSSCT00000033825 transcript, but not in ENSSSCT 
00000017641 (Table 2). Thus the promoter SNPs may 
affect the expression of a specific transcript, and together 
with amino acid changes, may influence the function of 
PRKAG3. Interestingly, the haplotype with the lowest 
P-value and a positive association with meat quality in both 
of the studied breeds is similar to a wild boar haplotype, 
with the exception that the wild boars' 199V is replaced 
by 1991 

Conclusions 

A single mutation in ASGA00070625, in the promoter 
region, in the amino acid at positions 24 or at 199 of 
PRKAG3 is not alone sufficient to create a favorable ef- 
fect on meat quality. Instead, a combination of variations 
or a haplotype with the promoter region variants g- 



157C and g.S8A and amino acid substitutions 24E, and 
1991 of PRKAG3 is needed to achieve a positive impact 
on meat quality traits, at least in the Finnish Yorkshire 
and Landrace populations. The results presented here 
can be directly applied in marker-assisted selection to 
improve the quality of pork meat. 

Methods 

Animal material for this study included previously col- 
lected semen and hair samples of the boars thus no eth- 
ical approval was required. All phenotypic data were 
kindly supplied by the Figen Ltd (http://www.figen.fi). 

Animals and meat quality measurements 

The study included 220 Finnish Yorkshire and 230 
Finnish Landrace AI (artificial insemination) boars. Add- 
itionally, three European wild boars were sequenced, but 
no phenotypic observations were available for these 
boars. Breeding values of the studied boars were esti- 
mated using the full national pig registry data including 
meat quality measurements from several thousand ani- 
mals. We used a single-trait BLUP procedure to estimate 
a breeding value for meat pH, color L* (lightness of 
meat) and a* (redness of meat). The statistical model in- 
cluded slaughter batch and sex as fixed, and litter and 
animal as random effects. The model was the same as 
used in national breeding value estimation in Finland, 
except that in the national evaluation all meat quality 
traits are analyzed simultaneously by a multitrait model, 
whereas in this study each trait was analyzed separately. 
We selected the single-trait approach to ensure that gen- 
etic correlation between traits did not affect the associ- 
ation results. The estimated breeding value (EBV) 
reflects the relative genetic merit of an animal. EBVs are 
generally more reliable than the animals own phenotype, 
because they are based on all available records on rela- 
tives and are simultaneously corrected for specific sys- 
tematic and non-systematic effects specified in the 
estimation model. Most of the meat quality data for a 
specific AI boar is obtained from its progeny and its full- 
and half-sibs. 

EBVs for meat quality traits are based on measure- 
ments taken from animals raised in a test station. Young 
piglets (on average 30 kg weight) are raised up to 100 kg 
live weight in a standardized test station environment. 
After the test period, all but the best young boars are 
sent to a slaughterhouse where pH and color L* and a* 
of meat are measured 24 h after slaughter. For this 
study, color L* and a* were measured on a freshly cut 
muscle surface with a Minolta CR 300 colorimeter and a 
CIELAB color scale standard [19,20], and pH was deter- 
mined using a Knick 752 pH meter and an Ingold 406 
electrode. Measurements were taken from loin (longissi- 
mus) and ham (semimembranosus) muscles. For more 
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information on the measurement procedures, see Sevon- 
Aimonen et al [1]. 

The studied Finnish Yorkshire and Landrace boars were 
born between 1992 and 2009, and included several relative 
pairs such as sire-son, full-sibs, grandsire-grandson, etc. 
Average relatedness between boars was 0.16 and 0.14 for 
the Yorkshire and Landrace, respectively. 

Genotyping and sequencing 

Part of the SNP data presented in this study originate from 
our previous whole-genome analyses [11,14] using the Por- 
cineSNP60 BeadChip (Illumina Ltd, San Diego, USA). 
Genotyping was performed at FIMM (Institute for Molecu- 
lar Medicine Finland, Helsinki, Finland) or at GeneSeek 
(Lincoln, USA). DNA was extracted either from hair folli- 
cles or semen, with a target DNA concentration of 300 ng. 
SNPs were mapped to the pig genome build Sscrofal0.2. 
We restricted our statistical analysis to cover only SNPs lo- 
cated in a 20-Mb region surrounding PRKAG3 (from 
120 Mb to 140 Mb on chromosome 15), because our previ- 
ous analyses had shown that the most significant SNPs for 
meat quality were in this region. 

After designing primer pairs for genomic sequence ana- 
lysis, we amplified the DNA fragments with gene-specific 
primers. PCR amplicons were purified using ExoSAP-IT™ 
(GE Healthcare, Piscataway, USA), and sequenced in both 
directions with the same primers as in the amplification 
procedures. Sequencing was performed on a 3500 x L Gen- 
etic Analyzer (Applied Biosystems, Carlsbad, USA) using a 
BigDye Terminator v3.1 kit (Applied Biosystems, Carlsbad, 
USA) and EtOH precipitation. 

Statistical method 

Prior to association analysis, the EBVs were deregressed 
and their weights were calculated by the method pro- 
posed by Garrick et al. [21]. The method removes parent 
average effects on EBV, so that the deregressed EBV 
more closely reflects the animals own performance and 
the performance of its offspring. The deregression pro- 
cedure also prevents regression towards the population 
mean, which is typical for EBVs which are based on a 
limited amount of information. Generally, the more reli- 
able the deregressed EBV is, the more weight it receives 
in the association analysis. 

Association analysis was performed either for individ- 
ual SNPs or a combination of SNPs (haplotype). Each 
SNP/haplotype was analyzed separately for association 
with meat quality traits using the following mixed linear 
model: 

y.=p + b*Xi + at + e h 

where y t is the deregressed EBV of the meat quality trait; 
%i is the number of minor alleles (0, 1, or 2) of the tested 



SNP or the number of copies of the tested haplotype (0: 
an animal carries no copies; 1: an animal carries one 
copy; 2: an animal carries two copies); b is the corre- 
sponding regression coefficient; a t is a random polygenic 
effect with a normal distribution with mean 0 and a 
variance-covariance structure of Acr 2 ^, where A is the 
additive relationship matrix and cr 2 ^ is the polygenic vari- 
ance; and e t is a random residual effect with a normal 
distribution with mean 0 and a variance-covariance 
structure of Io^/Wj , where I is an identity matrix, (? e is 
the residual variance, and w t is the weight. Association 
analyses were performed using the AI-REML method in 
the DMU program package [22]. Haplotypes were esti- 
mated with FastPHASE [23], and linkage disequilibrium 
plots were produced with Haploview [24]. 
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